#include "Acomb_func_with_acut.h"
bool Combine2_direct(TLorentzVector mom1, TLorentzVector mom2, ACuts &cuts, float &re_value){

	double mass_comb=TLV_M(mom1, mom2);
	double mt_comb  =TLV_Mt(mom1, mom2);
	double pt_comb  =TLV_Pt(mom1, mom2);
	double eta_comb =TLV_Eta(mom1, mom2);

	// mass window cut
	if(cuts.m.SWITCH){
		if( mass_comb< cuts.m.d  || mass_comb > cuts.m.u){
			return(false); 
		}
	}
	if(cuts.mt.SWITCH){
		if( mt_comb< cuts.mt.d  || mt_comb > cuts.mt.u){
			return(false); 
		}
	}
	if(cuts.pt.SWITCH){
		if( pt_comb< cuts.pt.d  || pt_comb > cuts.pt.u){
			return(false); 
		}
	}
	if(cuts.eta.SWITCH){
		if( eta_comb< cuts.eta.d  || eta_comb > cuts.eta.u){
			return(false); 
		}
	}
	re_value=mass_comb;
	return(true);
}
